data <- read.table("http://www.statpower.net/Content/312/Lecture Slides/manova.txt",header=T)
attach(data)
data
Y <- cbind(tear,gloss,opacity)
fit <- manova(Y ~ rate * additive)
### Multivariate Tests
summary(fit, test="Roy")
### Univariate Tests
summary.aov(fit)